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Abstract. South-east Asian mammals face a particularly severe threat of extinction. Borneo, the third largest island 
in the world, is located in the centre of South-east Asia. It harbours more endemic carnivores than does any other 
island except Madagascar. Almost half the Bornean carnivore species have been classified by The IUCN Red List 
of Threatened Species as threatened. Because little is known about most Bornean carnivores, predicting their spatial 
distribution is important for management strategies to improve the conservation of these species. As a part of the 
1 st Borneo Carnivore Symposium (BCS) we started to assemble a knowledge base of Bornean carnivores. We 
established the Borneo Carnivore Database which contains the previously largely fragmented occurrence records of 
carnivores on the island and then used these records to predict the distribution of 20 Bornean carnivores (all native 
species except sun bear Helarctos malayanus and the four otter species, Eurasian otter Lutra lutra, Asian small- 
clawed otter Aonyx cinereus, hairy-nosed otter Lutra sumatrana and smooth-coated otter Lutrogale perspicillata). 
We describe general considerations - the underlying assumptions, advantages, and most importantly the limitations 
and constraints - of species distribution modelling. We then summarise the methodological framework of our 
modelling approach and results of the sensitivity analyses. We emphasise that despite the extensive efforts to 
compile existing information, so few or spatially biased occurrence records exist for some species that the model 
outcomes presented in this journal issue must be interpreted cautiously. We recommend using new data as they 
become available to test our projections and improve our understanding of carnivore distributions on Borneo. 
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Abstrak (Bahasa Indonesia). Mamalia Asia Tenggara menghadapi ancaman dari kepunahan yang demikian tinggi. 
Borneo, sebagai pulau ke tiga terbesar di dunia yang terletak di tengah wilayah Asia Tenggara, menjadi tempat bagi 
banyak kamivora endemik dibandingkan pulau lainnya, terkecuali Madagaskar. Hampir setengah dari kamivora 
Borneo masuk klasifikasi Terancam dalam daftar IUCN Red List of Threatened Species. Mengingat masih sangat 
terbatasnya pengetahuan mengenai kebanyakan jenis kamivora Borneo, pendugaan penyebaran wilayah menjadi 
penting dalam strategi pengelolaan untuk peningkatkan konservasi dari jenis-jenis yang ada. Sebagai bagian dari 
hasil Simposium Kamivora Borneo Pertama (the 1st Borneo Carnivore Symposium; BCS) kami mulai menyusun 
pengetahuan dasar yang ada mengenai kamivora tersebut. Kami membangun Pangkalan Data Kamivora Borneo 
yang memuat sebagain besar infomiasi peijumpaan yang terpencar mengenai kamivora di Borneo dan menggunakan 
catatan tersebut untuk memperkirakan penyebaran 20 kamivora Borneo (semua satwa asli setempat terkecuali 
Beruang Madu Helarctos malayanus dan empat jenis berang-berang/sero yaitu Berang-Berang Pantai Lutra lutra, 
Sero Ambrang Aonyx cinereus, Berang-Berang Sumatera Lutra sumatrana serta Berang-Berang Wregul Lutrogale 
perspicillata. Kami menjelaskan pertimbangan umum - asumsi yang menjadi dasar, keuntungan, dan terpenting 
keterbatasan dan tantangan - pemodelan persebaran jenis tersebut. Kami kemudian menyimpulkan kerangka metode 
dari pendekatan pemodelan kami dan hasil dari analisa kesensitifitasan. Kami tekankan, walaupun usaha besar untuk 
mengumpulkan informasi yang ada telah berjalan, ada sedikit atau terdapat data spasial catatan keberadaan yang 
bias pada beberapa jenis, sehingga model yang disajikan dalam jumal ini hams diinterpretasikan dengan kehati- 
hatian. Kami merekomendasikan penggunaan data baru, ketika sudah tersedia, untuk menguji kembali proyeksi 
kami dan meningkatkan pemahaman kita mengenai persebaran kamivora di Borneo. 

Abstrak (Bahasa Malysia). Mamalia Asia Tenggara menghadapi ancaman kepupusan yang begitu tinggi. Borneo, 
sebagai pulau ketiga besar di dunia dan terletak di pertengahan Asia Tenggara, adalah tempat di mana terdapat paling 
banyak kamivora endemik (tidak didapati di tempat lain) dibandingkan kepulauan lain di dunia kecuali Madagascar. 
Hampir separuh daripada kamivora ini tersenarai sebagai terancam di atas IUCN Red List of Threatened Species. 
Oleh kerana pengetahuan kami sangat terhad mengenai kebanyakan kamivora Borneo, ramalan tentang penyebaran 
spesis-spesis ini adalah penting untuk pembentukan strategi pengelolaan dan meningkatkan usaha pemuliharaan 
mereka. Sebagai satu hasil dari Simposium Kamivora Borneo Pertama (1st Borneo Carnivore Symposium), kami 
mula berkumpul pengetahuan asas tentang kamivora ini. Kami mewujudkan satu Pangkalan Data Kamivora Borneo 
(Borneo Carnivore Database) untuk berkumpul maklumat tentang kamivora Borneo dan meramalkan penyebaran 
bagi 20 spesis kamivora Borneo (semua spesis yang didapati di Borneo kecuali Beruang Madu Helarctos malayanus 
dan empat spesis memerang iaitu Memerang ‘Eurasian’ Lutra lutra, Memerang Kecil Aonyx cinereus, Memerang 
Kumis Lutra sumatrana dan Memerang Licin Lutrogale perspicillata. Kami menjelaskan pertimbanganan am - 
penyangkaan dasar, faedahnya, dan yang paling penting, kebatasan dan masalah - menggunakan cara pemodelan 
penyebaran spesis jenis ini. Kami kemudian meringkaskan cara bagaimana pendekatan pemodelan ini dilakukan dan 
keputusan analisa kepekaan (sensitivity analysis). Kami menekankan bahawa walaupun usaha besar telah dilakukan 
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untuk mengumpulkan maklumat terkini, rekod-rekod yang sedia ada terlalu sedikit atau terdapat kecondongan 
setempat (spatial bias) untuk beberapa spesis sehingga keputusan model yang dikemukakan di dalam terbitan 
jurnal ini hams diinterpretasikan dengan hati-hati. Kami cadangkan agar bila data barn sudah ada, ianya digunakan 
untuk menguji peramalan model kami dan meningkatkan pemahaman mengenai penyebaran karnivora di Borneo. 


PROJECT BACKGROUND 

The 2008 Global Mammal Assessment of the threat status 
of the world’s mammal species by the International Union 
for Conservation of Nature (IUCN) revealed that South-east 
Asian mammals face a particularly severe threat of extinction 
(Schipper et al., 2008; Hoffmann et ah, 2011). On Borneo, 
almost 50% of carnivore species are listed as threatened on 
The IUCN Red List of Threatened Species. Because little 
is known about most Bornean carnivores, predicting their 
spatial distribution is important for developing management 
strategies to conserve these species (e.g., Kanagaraj et 
al., 2013). Species distribution models (SDM) relate 
environmental variables and species occurrence records to 
predict the distribution of species across large spatial scales 
(e.g., Elith & Leathwick, 2009). While spatial modelling 
studies stress the importance of predicting future species 
distributions in changing environments (e.g., Struebig et 
ah, 2015a, b), it is also important to establish historical 
distributions of species to provide baseline conditions for 
understanding recent shifts in distributions (Louzao et ah, 
2013). 

Our aim was to assist conservation efforts for Bornean 
carnivores by predicting potential historical and current 
distributions of carnivores (all native species except sun 
bear Helarctos malayanus (Raffles) and the four otter 
species Eurasian otter Lutra lutra (Linnaeus), Asian small- 
clawed otter Aonyx cinereus (Illiger), hairy-nosed otter 
Lutra sumatrana (Gray) and smooth-coated otter Lutrogale 
perspicillata (Geoffroy Saint-Hilaire)) throughout the island 
and the effects of human-dominated landscapes (habitat 
suitability). This is necessary to assess areas in which 
suitable habitat has potentially been lost. To this end, we 
first assembled the available species occurrences in a digital 
database. We then followed the hierarchical modelling 
approach by Wilting et ah (2010) where we estimated the 
potential historical distributions of 20 carnivore species 
based on environmental data and the suitability of current 
land area (land-cover class and human population density) 
using respondents’ knowledge and thoughts to calculate the 
current potential distribution based on a habitat suitability 
index (HSI) integrating historical distributions and current 
land-cover (Fig. 1). The modelling methods are also published 
in Kramer-Schadt et ah (2013) and Struebig et ah (2015a, b). 
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GENERAL CONSIDERATIONS REGARDING 
SPECIES DISTRIBUTION MODELS (SDMS) 

Maximum entropy modelling (program MaxEnt; Phillips et 
ah, 2006) has become popular for predicting distributions of 
rare and threatened species as it purportedly accommodates 
sparse data irregularly observed or sampled or low locational 
precision, typical for rare and elusive species, museum data, 
or species occurring in difficult terrain and climate (Elith et 
al,. 2006, 2011). Because systematic surveys for carnivores 
across Borneo have not been conducted, we used MaxEnt 
to estimate distributions that would account for spatially 
uneven data, i.e. data where the survey or search effort, 
respectively, was not equal everywhere. 

Despite the advantages of MaxEnt, there are several 
limitations of SDMs (Araujo & Guisan, 2006; Yackulic et 
ah, 2013). Researchers must be mindful of study design 
and biases inherent in presence data, and be cautious when 
interpreting model predictions (Pearce & Boyce, 2006). That 
is, MaxEnt is a presence-only approach (Li & Guo, 2011; 
Royle et ah, 2012) that because of methodological constraints 
(see Yackulic et ah, 2013) precludes the comparison of 
MaxEnt outputs across species. The output of MaxEnt 
models is a relative occurrence probability in relation to a 
default occupancy probability and cannot be compared across 
species observed during different temporal and spatial scales 
(Elith et ah, 2011), because occupancy probability differs 
across species. A second limitation using presence-only data 
is that sample selection bias, i.e., bias created by uneven 
search-efforts or non-systematic surveys, has a stronger 
effect on presence-only models than on presence-absence 
models (Phillips, 2008; Phillips & Dudik, 2008; Phillips et 
ah, 2009; Elith et ah, 2011). Sample selection bias occurs 
when outside a structured survey framework some areas in 
the landscape are surveyed more intensively than are others, 
as is common with museum data or opportunistic occurrence 
records such as sightings or road kills (Graham et ah, 2008). 
This uneven search-effort is often termed ‘sampling bias’. 

Our occurrence dataset was strongly biased to a few areas in 
Sabah, north-eastern Borneo, making the modelled relative 
probability of occurrence partly dependent on searching or 
surveying intensity (see Kramer-Schadt et ah (2013) for a 
detailed discussion on the effects of sampling bias). Although 
we reduced search-effort bias by filtering the occurrence 
records spatially and adjusting the background maps (see 
below and Kramer-Schadt et ah, 2013; Struebig et ah, 
2015a), this bias still affected model results. Consequently, 
we describe how results in the following papers of this 
supplementary issue should be interpreted. Models derived 
from data with large search-effort biases (e.g., few and 
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Fig. 1. Overview of the hierarchical modelling approach to develop distribution models for Bornean carnivores. In step one, the species 
occurrence records, environmental input layers and bias files representing record densities were used in a MaxEnt modelling approach to 
derive historical habitat suitability or possible bioclimatic distribution. In step two, the historical habitat suitability map was integrated 
with reclassified land-cover and human population density suitability scores (habitat associations) in a habitat suitability index model to 
yield relative current habitat suitability probability or potential current distribution. 


spatially clumped records for some species in Sabah; see Fig. 
2), generally produce reliable predictions for areas intensively 
surveyed but not for other areas. Large areas with few 
records (e.g., parts of West Kalimantan) should be regarded 
as unstudied and areas with low predicted probability of 
occurrence in these areas should not be considered necessarily 
to be unsuitable. Underrepresentation of non-surveyed areas 
in the respective species’ distribution maps can further be 
reinforced when the climatic conditions are very much 
different from areas from whence locations stem. Hence, areas 
in Brunei, being similar to climatic conditions in Sarawak 
and Sabah, might be predicted correctly, whereas data-poor 
areas that are climatically different from areas with records, 
e.g., in South Kalimantan (Fig. 3A), cannot be predicted 
as suitable. Thus, care must be taken when discussing low 
suitability in the predictions as to whether this is because of 
methodological issues, climatic differences or because large 
areas have been transformed into unsuitable land-cover, e.g., 
oil palm plantations (Fig. 3B). Here, a careful evaluation 
of the model results in the context of survey effort and the 
individual species’s biology is needed. 


Also, because of search-effort bias the comparison of 
absolute model prediction values across regions for the same 
species is not meaningful. That is, a relative probability of 
occurrence of 0.4 in Kalimantan (few records) could be 
equivalent to value of 0.7 in Sabah (many records). Finally, 
because each model is based on the records of the respective 
species, model predictions are based on a species-specific 
scaling. However, the detection probability is fixed at the 
same threshold for each species. Thus, comparisons across 
species are not possible, i.e. species A is not necessarily 
more likely to occur within a certain area than does another 
species B even if species A’s modelled relative occurrence 
probability in that area is higher than that for the species B. 

Despite these limitations, the SDMs provide an important 
contribution towards knowledge of Bornean carnivores. This 
is the first time that distribution records about all Bornean 
carnivores (except sun bear and otters) were assembled and 
processed (i.e., modelled) in a standardised framework; 
this provides a baseline for compiling additional data. The 
SDMs provide an explicit quantification of habitat suitability 
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Fig. 2. Location of Borneo carnivore occurrence records by record-type for spatial precision categories 1-4 including countries as well 
as state boundaries. 


assessments, can be refined with new records, and can 
serve to facilitate discussion of species conservation in a 
spatial context. Further, they can be used to identify areas 
with limited information for further research as well as 
conservation areas and potential corridors. Finally, SDMs 
can function as a tool to test ecological hypotheses about 
species’ distributions and niches and to help synthesise 
existing knowledge. 

METHODOLOGICAL APPROACH 

I: THE BORNEO CARNIVORE DATABASE 

The Borneo Carnivore Database was created to assemble 
occurrence records for the 20 Bornean carnivores, the subject 


of this proceedings issue (see also Table 2). All records with 
uncertain species identification or with uncertain or missing 
geographic information were excluded from the database. We 
collected and assembled occurrence records from scientists, 
wildlife observers, museums and literature (Table 1). We 
also reviewed the Indo-Australian Mammal Atlas (van Strien, 
2001) for occurrence records. Data types included museum 
specimens found in the Mammal Networked Information 
System (MaNIS) and Global Biodiversity Information 
Facility (GBIF) non-fossil museum databases or provided by 
the museums directly, road kills, hunters’ relicts like skins 
and skulls where the spatial origin was established, camera¬ 
trapping data and direct sightings. Indirect information in the 
forms of tracks and other signs, and information provided 
by villagers and other local people was excluded because of 
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Fig. 3. A, climatic differences in Borneo expressed here as the Driest Quarter (Flijmans et al, 2005; www.worldclim.org/bioclim; bioclim 
17). B, land-cover map used for habitat associations. 
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the high potential of species misidentification (e.g., Mohd- 
Azlan et al., 2013). 

Each occurrence record was validated with other records and 
those in question were discussed with respective contributors. 
Contributors were asked to provide only records of high 
certainty and camera-trapping photographs were shared with 
a broader community for species identification as required. 
Due to our careful review of the records and our strict criteria, 
most uncorroborated records from public databases such as 
MaNIS, GB1F and the lndo-Australian Mammal Atlas were 
excluded from modelling because of missing or uncertain 
information, but kept as entry in the database. 

Each record included the following information when 
available: species, geographic coordinates in decimal degrees, 
location name, collector or observer, date, type of record, 
museum identification, museum name, habitat as provided 
by the observer, altitude as provided by the observer, 
reference and remarks. When geographic coordinates were not 
available, we assigned coordinates using maps of Kalimantan 
(CV. Indo Prima Sarana; http://www.indoprimasarana.com/), 
Google Earth and knowledge from researchers working in 
those regions. 

We categorised records based on spatial precision provided 
by the contributors of (1) within 0.5 km, (2) 0.5-2.0 km, 
(3) 2.0-5.0 km, (4) over 5.0 km and (5) unknown radius 
(Table 2, Fig. 2). Overall, we had compiled 3979 records 
by September 2011, of which 1151 were from 36 museums, 
2290 were from 46 field studies, 126 were from online 
databases, and 412 were found in 67 of 270 documents 
published from 1845 to 2011 (Tables 1, 2). Duplication of 
records was avoided by cross-validating the collector’s name 
or location specification. Flowever, if there were differences 
in the sources, i.e., entry in a museum list, book or database, 
we could not distinguish the record at this stage. Flowever, 
if this occurred, it would not have affected the modelling 
procedure because such records would have been spatially 
too imprecise or spatially filtered (see below). Of these, 1850 
records were gathered before 2001, 2078 records were from 
2001 or later, and 51 lacked the date. From field sources, 
741 records were sightings and 1615 were camera-trap 
photographs. From a total of 3581 records in categories 1—4 
we compiled 2464 (69%) records from Malaysian Borneo 
(578 [16%] from Sarawak and 1886 [53%] from Sabah), 
1100 (31%) from Indonesian Kalimantan (West Kalimantan 
177 [5%], South Kalimantan 29 [1%], Central Kalimantan 
336 [9%], East Kalimantan 389 [11%], North Kalimantan 
169 [5%]) and 17 (0.5%) from Brunei Darussalam. 

Species occurrence records used for modelling. We 

only used records with spatial precision categories 1-3 for 
generating SDMs because of potential spatial inaccuracy 
of records with precision categories 4 and 5. To reduce 
inaccuracy in model projections associated with spatial 
autocorrelation in model residuals because of a clustered, 
dependent data structure (cf. Dormann et al., 2007; Veloz, 
2009), we applied further spatial filtering and used only one 
record within a radius of 2 km (greater than the approximate 


home range radii of many small carnivores) for each species. 
For species with a greater number of records this radius 
was extended to 5 or 10 km (Table 3), to reduce possible 
search-effort bias without removing too many data points. 
We retained records with greatest location accuracy; when 
data clusters included two or more records of equal accuracy, 
we selected the record to be retained randomly. This 
spatially filtered dataset (M 2 in Tables 1 of the respective 
species papers) was used to model only the distribution 
within the Malaysian state of Sabah (referred to as ‘spatial 
filtering model’ M 2 ). Flowever, even after spatial filtering, 
the Malaysian State of Sabah, which represents 10% of the 
area of Borneo, still contained over 50% of all records. To 
model the distribution outside Sabah appropriately, we further 
reduced the number of records for each species by selecting 
randomly the number of records in Sabah which corresponded 
by areal proportion to the number of records outside Sabah 
(balanced dataset M, in Tables 1 of the respective species 
papers in this supplementary issue, referred to as ‘balanced 
model’ M,). For example, if 90 records of a species were 
detected outside Sabah, we included 10 records of that 
species from Sabah from the respective species database 
(see Kramer-Schadt et al., 2013, for sensitivity analyses of 
spatial filtering effects). 

II: THE BORNEO GEODATABASE 

All analyses were done using Arclnfo 9.3 and ArcGIS v. 
9.3.1 (ArcView, Spatial Analyst, 3D Analyst, Geostatistical 
Analyst) by Environmental Systems Research Institute (ESRI) 
and R 2.10.1 (R Development Core Team 2009). All maps 
were geo-referenced, clipped or shifted to the same extent 
with a cell size of 30 arc seconds (~ 1 km 2 ). We filled 
missing data in the Bioclim maps (see below) by applying 
mean values of surrounding cells using a moving window 
approach across a radius of eight cells around the focal cell 
(procedure ‘focalmean’ in Arclnfo). 

Environmental input variables: maps for modelling the 
potential historical distribution. We used 19 global climate 
data maps including temperature and precipitation parameters 
calculated from long-term time series (1950-2000) of 
monthly values (Hijmans et al., 2005; http://www.worldclim. 
org/bioclim; Fig. 3A; for a detailed list see Table SMI 
subtext). We also created maps that contained information 
on landscape attributes such as aspect, terrain ruggedness 
and distance to water courses, which are important for 
some species such as flat-headed cat Prionailurus planiceps 
Vigors & Horsfield (Wilting et al., 2010). Aspect and terrain 
analyses were based on a digital elevation model (DEM, 
90m raster resolution; http://srtm.csi.cgiar.org). Aspect was 
calculated using the extension Spatial Analyst in ArcView 
9.3.1. The Topographic Ruggedness Index (TRI, Riley 
et al., 1999) expresses the elevational difference between 
each cell of the DEM and its eight nearest neighbours 
(Moore neighbourhood), with ruggedness classes ranging 
from one (level) to seven (extremely rugged). Hydrologic 
analyses (Hydrologic Analysis in Spatial Analyst) included 
flow directions calculated to derive a subsequent river net. 
We used a single flow direction approach, that is, the flow 
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Table 1. Reference numbers for museums and for individual contributors. There is no connection between the Museums (left) and the 
individual (right). 


Museum Name 


Scientists 


American Museum of Natural History, New York, U.S.A. (AMNH) 

Museum fur Naturkunde, Berlin, Germany (ZMB) 

Bernice P. Bishop Museum, Honolulu, Hawaii, U.S.A. (BPBM) 

Natural History Museum, London, U.K. (BMNH) 

Museum Zoologicum Bogoriense, Bogor, Java, Indonesia (MZB) 

Zoologisches Forschunginstitut und Museum Alexander Konig, Bonn, Germany (ZFMK) 

Brunei Forest Center, Sungai Liang, Brunei, Darussalam 

Brunei Museum, Bandar Seri Begawan, Brunei Darussalam 

Institut Royal des Sciences Naturelles de Belgique, Brussels, Belgium (IRSNB) 

Staatliches Museum fur Tierkunde, Dresden, Germany (MTKD) 

Field Museum of Natural History, Chicago, Illinois, U.S.A. (FMNH) 

Naturmuseum und Forschungsinstitut Senckenberg, Frankfurt, Germany (SMF) 
Zoological Survey of India, Calcutta, India (ZSI) 

University of Kansas Museum of Natural History, Lawrence, Kansas, U.S.A. (KU) 
University of Malaya, Kuala Lumpur, Malaysia (UMKL) 

Musee National d‘Histoire Naturelle, Paris, France (MNHN) 

Museo Civico di Storia Naturale Giacomo Doria, Genova, Italy (MSNG) 

Museum of Comparative Zoology, Harvard University, Cambridge, Massachusetts. 
U.S.A. (MCZ) 

University Museum of Zoology, University of Cambridge, U.K. (UMZC) 

National Museums of Scotland, Edinburgh, Scotland (NMSZ) 

Universitetets I Oslo, Zoologisk Museum, Oslo, Norway (ZMUO) 

Rijksmuseum van Natuurlijke Historie ‘Naturalis’, Leiden, Netherlands. (RMNH) 
Naturhistorisches Museum Basel, Basel, Switzerland (NMBA) 

Naturhistorisches Museum, Vienna, Austria (NMW) 

OZCAM (Online Zoological Collections of Australian Museums) Provider 
Sabah Museum, Kota Kinabalu, Sabah, Malaysia (SBM) 

Collection of Sabah Parks, Sabah, Malaysia 
Sarawak Museum, Kuching, Sarawak, Malaysia (SMK) 

National Museum of Natural History, Washington D.C., U.S.A. (USNM) 

Staatliches Museum filer Naturkunde in Stuttgart, Germany (SMNS) 

Naturhistoriska Riksmuseet, Goteborg, Sweden (NHMG) 

National Museum of Victoria, Melbourne, Victoria, Australia (NMV) 

Universiteit van Amsterdam, Zoologisch Museum, Amsterdam (ZMA) 

Zoological Museum, University of Copenhagen, Copenhagen, Denmark (ZMUC) 
Zoological Reference Collection, National University of Singapore, Singapore. (ZRC) 
Zoologische Museum Staatssammlung, Miinchen, Germany (ZSM) 
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Azlan Mohamed & Andreas Wilting* 

Ben Wielstra & Sander Pieterse* 

Brent Loken 
Craig Robson 
Danielle Kreb 
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Edward Pollard 
Erik Meijaard 
Gabriella Fredriksson 
Godwin Limberg 

Gono Semiadi & Hiromitsu Samejima* 
Hans Skotte Moeller 
Henry Bernard 
Hiromitsu Samejima 
Hiroshi Sasaki 

J. W. Duckworth 
James Eaton 
Jason Hon 

Jedediah Brodie & Anthony Giordano* 
Joanna Ross & Andrew Hearn* 

John D. Pilgrim 
John Howes 
John Mathai 
Jon Hah 
Konstans Wells 
Lim Boo Liat 
Mark Bezuijen 
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Matt Heydon 
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Norman Lim 
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Table 2. Borneo carnivore occurrence records with spatial precision categories 1-3 (< 5 km radius) collected before 2001 and during 2001-2011; for 29 records no date could be assigned 
(abbreviations: ct = camera-trapping, ds = direct sighting, hr = hunters’ relic, rk = road kill). Literature included the van Strien (2001) database; scientist records include verified records collected 
in local communities and villages. Individual species localities provided by multiple sources were treated as described in the text. 
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Herpestidae 

Herpestes 

brachyurus 

Herpestes 

semitorquatus 





Table 2. Continued 
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from each cell was directed only to the steepest down-slope 
neighbouring cell. The cells in which flow accumulated were 
considered rivers according to the number of accumulated 
cells. We established flow accumulation cut-offs to 5000, 
100,000 and 500,000 cells for small, medium and large 
perennial river systems, respectively. Corresponding 
catchment sizes were about 42 km 2 , 850 km 2 and 4230 km 2 , 
respectively. By combining the three river net maps with 
a map of the water bodies of Borneo (i.e., lakes and large 
rivers; country inland water data of the Digital Chart of the 
World, downloaded from: http://www.diva-gis.org) three 
distance-to-water maps were created. 

Accounting for search-effort bias. Spatial distributions of 
records across Borneo were heterogeneous, therefore, we 
developed species-specific ‘bias files’ representing relative 
search efforts. We first mapped species records on our 
raster-based grid, and then assigned each cell with a score 
of 1 for species presence or 0 for absence. We summed the 
values within each cell’s Moore neighbourhood to generate 
a map of search-effort density. We assigned the value 0.1 to 
cells with scores of 0, representing a 10-fold lower search 
effort than cells with a value of 1 (Kramer-Schadt et al. 
(2013) gave a detailed sensitivity analysis regarding bias 
file background values). 

Habitat associations: maps for modelling the potential 
recent distribution. For each species, the suitability of 
land-cover classes and human population density to serve 
as each species’s habitat was classified from 0 to 4 based on 
respondents’ opinions, then rescaled from 0 to 1 (see below 
and Tables 2 in respective species papers). The respondents 
were scientists and conservationists working on carnivores 
and other species on Borneo. We calculated mean suitability 
score values from the respondents’ scoring per map and 
species. From these averaged respondents’ rankings we 
reclassified the land-cover attributes for species-specific (i) 
land use (Z,) and human population density (P,) suitability 
maps. 

Human population density (LandScan 2007™ High 
Resolution Global Population Data Set, Oak Ridge National 
Laboratory, UT-Battelle, LLC) was categorised as 0 = 
population densities too high for the species to occur regularly 
and therefore unsuitable, 0.25 = large negative effects, 0.5 
= medium effects, 0.75 = little impact and 1 = no effect 
on the species. Then, the mean values were calculated 
using the respondents’ rankings. Despite some differences 
between respondents in ranking human population density 
as suitable habitat for each species, the influence of this 
variable on model results was quite low because high human 
population density was strongly associated with modified 
land-cover classes such as plantations, mixed crops, degraded 
forests, bare areas and burnt areas. Thus, information about 
respondents’ reclassifications of human population density 
is omitted from the single-species papers, which focus on 
land-cover reclassification instead. 
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For land-cover, we used the PALSAR land-cover map 
validated for Borneo (Hoekman et al., 2009) from 2007 
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Neofelis diardi 219 94 5 48 78 0.772 0.569 

Pwdofelis marmorata 80 69 2 27 66 0.838 0.589 

Catopuma badia 47 40 2 23 51 0.775 0.687 

Prionailurus planiceps 105 76 2 50 52 0.789 0.713 

Prionailurus bengalensis 297 102 10 63 79 0.701 0.645 
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Table 4. Reclassification of the PALSAR land-cover map incorporating a digital elevation model (DEM) (see also Fig. 3 B). 


Original Classification 

Amended Classification 

Tropical lowland forest 

Rowland forest 0-500 m 

Upland forest 501-1000 m 

Tropical mountain forest 

Lower montane forest 1001-1500 m 

Upper montane forest above 1500 m 

Forest mosaics / fragmented or degraded forests 

Lowland forest: 

Forest mosaics/fragmented or degraded forests 

0-500 m 

Upland forest: 

Forest mosaics/fragmented or degraded forests 

501-1000 m 

*Lower montane forest: 

Forest mosaics / fragmented or degraded forests; 

1001-1500 m 
*Upper montane forest: 

Forest mosaics / fragmented or degraded forests; 
above 1500 m 

Riverine forest 

Swamp forest 

Swamp forest 


Peat swamp (pole) forest 


Peat swamp (padang) forest 


Mangrove forest 

Mangrove 

Nipah mangrove forest 


Plantations 

Old plantations 

High shrub 

Young plantations and agriculture 

Medium shrub 


Ferns/grass 

Burned forest area 

Grasslands 


Tree cover, burnt 


Cropland (upland) 

Mixed crops 

Cropland (irrigated) 

Water and fishponds * ** 

Water bodies 

Water bodies 

No data 

No data 


*Categories not applicable to the reclassified map. 

**Groundtruthing revealed irrigated cropland to be mainly fishponds or other human-modified wetlands and was therefore reclassed to 
‘water and fishponds’. 
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with an original resolution of 50 m and resampled to the 
cell-size of the bioclimatic predictors using the ‘majority’ 
approach in ArcGIS 9.3.1. We modified it with height 
information from the DEM in 500 m steps to allow a finer 
grain of land-cover classification (Table 4, Fig. 3 B). The 
respondents then assigned one of five categories to each 
land-cover type ranging from 0 to 1 with 0 = unsuitable 
habitat, 0.25 = low suitability, 0.5 = medium suitability (e.g., 
potential use as corridor), 0.75 = high suitability and 1 = 
very high suitability. For certain species some respondents 
ranked particular land-cover classes as unsuitable whereas 
others ranked them as highly suitable. This indicates how 
little is known for some Bornean carnivores. Table 2 in 
each species paper of this supplementary issue gives the 
scoring of the land-cover classifications and including the 
range of the ranking as an indicator of consistency among 
the respondents. 

Ill: THE BORNEO CARNIVORE DISTRIBUTION 
MODELLING 

Modelling the potential bioclimatic distribution. We 

used MaxEnt v. 3.3.3a (www.cs.princeton.edu/~schapire/ 
maxent; Elith et ah, 2006; Phillips et ah, 2006) to associate 
occurrence records with abiotic environmental conditions. 
We ran MaxEnt with the default settings and evaluated the 
percentage contribution of each variable with a jackknife 
procedure. We ran 10 replicates and used mean probabilities 
predicted for further analyses. As a measure of SDM accuracy 
or discriminative power, respectively, generally the threshold 
independent and prevalence insensitive presence-only 
(PO) area under the curve (AUC P0 ) produced by MaxEnt 
classifying background versus presence points are used 
(Yackulic et ah, 2013). Models with an AUC > 0.7 are said 
to have fair discriminative power (Hosmer & Lemeshow, 
1989); however, this is a statistical measure of how well the 
model could discriminate spatial variables depending on the 
input data. AUC and AUC PO are not always good predictors 
of model quality in terms of ecological plausibility (Lobo et 
ah, 2007; Kramer-Schadt et ah, 2013; Yackulic et ah, 2013). 
Therefore AUC PO are presented only in Table 3 and are used 
for our sensitivity analysis (see below) to compare model 
performance, but are not used as a measure of model quality. 

Habitat Suitability Index (HSI): modelling potential 
current distributions. Using MaxEnt, we produced a map 
estimating the relative probability of occurrence of each 
species i per grid cell (M.) that was updated using the current 
species-specific land-cover (Z,) and population density (P) 
maps to derive a habitat suitability index (HSI,) for potential 
recent distributions (Fig. 1) following a modified equation 
of Wilting et ah (2010) that places greater weight on the 
current land-cover map: 

HSI, = (M 2 * Z, 3 * Pi ) 1/6 (1). 

We placed lesser weight on the human population density 
map as species’ responses to human population densities 
are largely unknown for Bornean carnivores. However, 
by including the human population density map with a 


lower weight in analyses we attempted to account for this 
potential effect. 

Sensitivity analysis of the MaxEnt model. We conducted 
uncertainty analyses on the MaxEnt potential bioclimatic 
distribution layer to account for potential effects of 
clumping of species records because of search-effort bias, 
background adjustment with bias files and multicollinearity 
of environmental input layers. To this end, we compared 
model runs conducted with the full set of input records 
(termed ‘no reduction’) with systematically filtered species 
records (i.e. 2-10 km reduction termed spatial filtering and 
reduction proportionally to the size of Sabah termed balanced 
design). Additionally, we assessed the sensitivity of model 
outputs by omitting the bias files and we ran all models with 
a reduced set of uncorrelated (Pearson’s moment correlation 
|/j < 0.75) environmental input layers (see Table SM2 subtext 
for the reduced set). 

We evaluated sensitivity based on changes in the AUC P0 
values and applied two cut-off values based on the relative 
probability of occurrence (>0.5 and >0.7) and calculated 
the total area above these thresholds in relation to the total 
area of Borneo. These cut-off values were used only to 
assess relative changes in spatial model predictions within 
the same data set, and not used as thresholds to categorise 
suitable and unsuitable areas (Tables SMI, SM2; see also 
Kramer-Schadt et al., 2013). 

Results of the sensitivity analysis. Multicollinearity of 
input layers did not strongly affect model outputs, because 
the predicted areas increased only marginally (Table SM2; 
see also Kramer-Schadt et al., 2013). The use of a bias file 
did not affect AUC PO s (Table SMI), though the extent of 
spatial predictions increased with the use of a background 
bias file in sensitivity scenarios (Table SM2). This increase 
was negligible in the balanced scenarios and core areas (areas 
above cut-off with relative occurrence probability >0.7). 
Reducing clumping of records using spatial filtering also 
did not affect model AUC PO s for most species compared 
with model results using the full data set (Table SMI). By 
contrast, balancing the number of records in Sabah relative 
to the total survey effort in Borneo reduced AUC PO s for 
most species up to one order of magnitude (Table 3). 
However, Kramer-Schadt et al. (2013) demonstrated that 
model quality, despite the reduced AUC P0 in the balanced 
model, was actually higher in these models because they 
showed lower omission and commission errors than did 
models based on the full dataset. For well-surveyed areas 
such as Sabah the models with spatial filtering showed the 
lowest number of omission and commission errors (Kramer- 
Schadt et al., 2013). Therefore, we divided our analyses. 
Fig. 2 in most of the specific species papers (excepting 
those species recorded solely or very largely from Sabah) 
shows an insert (Figure 2A) for Sabah, which is based on 
the spatial filtering model, whereas Figure 2B is based on 
the balanced model. Each species distribution map in the 
respective species papers shows the protected area network, 
consisting of fully protected areas as well as, for Indonesia, 
forest reserves (such information was not available for 
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Brunei and Malaysia), because forests still can contribute to 
conservation if well managed (http://gislab.cifor.cgiar.org/ 
wm/borneo/; Struebig et al., 2015a). We updated this map 
with latest information on recently protected areas in Sabah, 
Malaysia, provided by the Sabah Forestry Department, Forest 
Resource Management (FRM) Division. 

Threshold maps for delineating possible conservation 
hotspots. To be able to compare all species models with 
their relative probabilities of occurrence (HSI value), we 
standardised all maps by converting the HSI values into 
binary maps separating the predictions into suitable and 
unsuitable. We used the commonly applied 10% and 25% 
omission error thresholds to yield a liberal and a conservative 
suitability prediction, respectively, allowing for some leeway 
in interpreting model results. These binary threshold maps 
were then combined in Mathai et al. (2016) to map those 
areas with the highest possible occurrence probability of 
selected (as indicated by the respective map) carnivore 
species. These combined overlays were used to discuss and 
delineate areas of conservation priorities. 

SYNTHESIS AND DISCUSSION 

To map predictions of the potential bioclimatic distributions 
(MaxEnt projections), we used the results of the balanced 
design M, (for Sabah, the spatial filtering design M 2 ) together 
with the bias file. The uncertainty analyses revealed that 
model AUCp 0 s were robust to multicollinearity and the use 
of bias files, however, AUC PO s were reduced below 0.7 for 
some species. This does not mean that models do not project 
species distributions effectively, because a perfect model can 
have a low AUC PO value if the species’s detection probability 
is low in a certain area or if it experiences frequent local 
extinctions (Yackulic et al., 2013). Although we considered 
potential sources of bias in our modelling approach and 
balanced our design as much as possible, spatial search- 
effort bias still affected model performance. That is, relative 
probability of occurrence produced by MaxEnt is partly 
related to search effort (see Kramer-Schadt et al., 2013). 
However, it is not possible with available data to account 
further for spatial bias. Consequently, species distribution 
maps should be interpreted with caution. 

To aid model validation and interpretation, a threshold is 
usually applied to develop binary maps of suitable and 
unsuitable areas (Liu et al., 2005). Severe variation in 
predicted suitable areas using varying thresholds can occur; 
thus the selection of a threshold can have a strong effect 
on results (Liu et al., 2005). To avoid this shortcoming, we 
used continuous relative probability measures to estimate 
species habitat suitability and used fixed thresholds only for 
sensitivity analyses (Table SM2). For the summarising map 
overlays (as in Mathai et al. 2016) we used the commonly 
applied thresholds (see above). 

The MaxEnt distribution was one of three factors we used 
to predict potential species distributions. A strength of our 
modelling approach is that we considered multiple sources 
of information to estimate habitat suitability. We placed 


greatest emphasis on current land-cover class that was refined 
by respondents’ knowledge. In doing so, we ameliorated 
model accuracy by updating predictions with rule-based 
models (Yamada et al., 2003). In line with our approach. 
Beck et al. (2014) showed that spatial filtering provided 
more accurate projections. 

Discussions regarding projected distributions, habitat 
associations and conservation needs of each species follow 
in the species-specific papers, where modelled relationships 
are examined using our current state of knowledge. That 
speaking, the maps do not predict the actual distribution of 
the species, they predict conditions suitable for the species 
in respect of the modelled attributes. Other non-modelled 
aspects like dispersal barriers and targeted overhunting may 
mean the species is not there. These aspects are thoroughly 
discussed in the single species papers. Our intent is that 
results of this effort will assist with future prioritisation of 
conservation and research efforts for carnivores on Borneo. 
New data, particularly additional occurrence records, can 
be used to refine the projected species distributions and 
provide independent assessments of projections presented 
in this initial effort. 
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SUPPLEMENTARY MATERIAL 

Table SMI. Sensitivity analysis of species distribution models of Bornean carnivores showing the three most important bioclimatic 
layers (Variable) with the highest contribution (% cont.) to explaining the distribution of records as well as AUC po . Only model results 
calculated from all listed environmental layers are shown, i.e. sensitivity analysis was not conducted for scenarios with ‘non-correlated’ 
environmental layers only. Bold: model runs (Balance models M,) used to assess current distributions in the habitat suitability index 
(HSI) throughout Borneo (except for Sabah projections which derive from the spatially filtered models M 2 ). For Bornean ferret badger 
Melogale everetti and Sunda stink-badger Mydaus javanensis only the projections of the Spatial Filtering M 2 model were used because 
of the range restriction of these species). See table subtext for detailed description of variables. 


Model 



Reduction Sabah (balanced) M, 



Bias File Used 


Yes 



No 


Species 

AUCpo 

Variable 1 

% Cont. 

AUCpo 

Variable 

% cont. 

Martes flavigiila 

0.67 

bio_02 

37.8 

0.671 

bio_02 

33.1 



bio_22 

17.4 


bio 22 

16.7 



bio_25 

12.9 


bio_25 

13 

Mustela nudipes 

0.763 

bio 07 

17.6 

0.755 

bio_04 

15.8 



bio 17 

12.8 


bio_07 

15.1 



bio_04 

12.2 


bio_19 

12.6 

Melogale everetti 

0.995 

bio_09 

28.7 

0.995 

bio_09 

24.8 



bio 12 

23.5 


bio 12 

24.1 



bio_17 

8.2 


bio_17 

9.2 

Mydaus javanensis 

0.746 

bio 02 

46.7 

0.751 

bio_02 

44.8 



bio_27 

34.3 


bio_27 

34.4 



bio_25 

6.1 


bio_25 

6.7 

Prionodon linsang 

0.65 

bio_18 

20.6 

0.648 

bio_18 

18.8 



bio_26 

17.2 


bio_26 

15.8 



bio_14 

9.3 


bio 14 

8.9 

Viverra tangalunga 

0.661 

bio_22 

30.7 

0.662 

bio 22 

29.7 



bio_08 

9.4 


bio_02 

9.1 



bio_02 

8.5 


bio_08 

8 

Paradoxurus hermaphroditus 

0.64 

bio_13 

14.2 

0.642 

bio_21 

13.3 



bio_21 

12.8 


bio_13 

12.9 



bio_15 

11.4 


bio 12 

12.1 

Paguma laivata 

0.69 

bio_15 

27.7 

0.686 

bio_15 

29.2 



bio_22 

20.8 


bio_22 

19 



bio_07 

8.5 


bio_13 

12.5 

Arctictis binturong 

0.701 

bio_07 

38.8 

0.707 

bio_07 

34.8 



bio_15 

17.5 


bio_15 

17.8 



bio_22 

12.9 


bio 22 

11.1 

Arctogalidia trivirgata 

0.645 

bio_02 

21.1 

0.637 

bio_14 

17.3 



bio_14 

19.5 


bio_13 

16.4 



bio_13 

13.1 


bio_02 

15.2 

Hemigalus derby anus 

0.642 

bio_07 

16.3 

0.643 

bio_13 

16.1 



bio_14 

15 


bio_07 

13.9 



bio_13 

13.6 


bio_14 

13.3 

Diplogale hosei 

0.957 

bio_05 

46.6 

0.856 

bio 22 

18 



bio_15 

26.5 


bio_09 

13.9 



bio_27 

8.1 


bio 12 

13.7 

Cynogale bennettii 

0.495 

bio_22 

19 

0.494 

bio 22 

18 



bio 19 

14 


bio_09 

13.9 



bio_09 

13.9 


bio_12 

13.7 
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Model 



Reduction Sabah (balanced) M, 



Bias File Used 


Yes 



No 


Species 

AUC P o 

Variable 1 

% Cont. 

AUCpo 

Variable 

% cont. 

Herpestes brachyurus 

0.668 

bio_15 

24.3 

0.666 

bio_15 

21.1 



bio 07 

16.2 


bio_07 

20.6 



bio 27 

12.8 


bio_27 

11.4 

Herpestes semitorquatus 

0.709 

bio_22 

30.7 

0.71 

bio 22 

29 



bio 27 

16.7 


bio_27 

19.3 



bio 07 

10.8 


bio_07 

8.5 

Neofelis diardi 

0.569 

bio_07 

25.3 

0.569 

bioJ)7 

23.6 



bio_26 

21.3 


bio 26 

20.5 



bio_22 

9.3 


bio 22 

8.8 

Pardofelis marmorata 

0.589 

bio 27 

15.4 

0.575 

bio_27 

16.2 



bio_26 

14.9 


bio 26 

14.2 



bio_02 

14.7 


bio_02 

13.8 

Catopuma badia 

0.687 

bio_22 

31.3 

0.686 

bio 22 

30 



bio 07 

21 


bio_07 

21.8 



bio_25 

10.6 


bio_25 

10.5 

Prionailurus planiceps 

0.713 

bio_13 

15.8 

0.716 

bio_13 

18.3 



bio 25 

15.2 


bio 25 

15.2 



bio_06 

14 


bio_09 

12.7 

Prionailurus bengalensis 

0.645 

bio_22 

23.6 

0.641 

bio 13 

22.3 



bio_13 

17.2 


bio 22 

19.9 



bio_02 

15.7 


bio_02 

16 

Sensitivity analysis - Alternative 1: 

no reduction of clumped data through spatial filtering 




Model 



No Reduction 



Bias File Used 


Yes 



No 


Species 

AUCpo 

Variable 

% Cont. 

AUCpo 

Variable 

% cont. 

Martes flavigula 

0.787 

bio_02 

40.3 

0.791 

bio 02 

27.9 



bio_04 

11.5 


bio_04 

11.7 



bio_08 

7.2 


bio_16 

8.7 

Mustela nudipes 

0.825 

bio 12 

17.2 

0.837 

bio_12 

17.1 



bio_02 

16 


bio_04 

12.2 



bio_04 

12.2 


bio 02 

11 

Melogale everetti 

0.996 

bio_09 

22 

0.997 

bio_12 

20.4 



bio 12 

17.6 


bio_09 

16.8 



bio_08 

13.4 


bio_08 

13.6 

Mydaus javanensis 

0.946 

bio_02 

37.3 

0.945 

bio 02 

29.3 



bio 18 

20.2 


bio_18 

17.4 



bio_04 

16.8 


bio 04 

16.2 

Prionodon linsang 

0.786 

bio 22 

27.4 

0.788 

bio 22 

23.9 



bio_18 

23 


bio_18 

18.3 



bio_04 

9.1 


bio 04 

11 

Viverra tangalunga 

0.86 

bio_02 

36.1 

0.862 

bio 02 

24.4 



bio_08 

12.1 


bio_08 

10.5 



bio_18 

11.3 


bio_18 

9.7 

Paradoxurus hermaphroditus 

0.839 

bio_02 

36.1 

0.846 

bio 02 

29.3 



bio 18 

17.2 


bio_18 

14.9 



bio_04 

15.8 


bio_04 

12.9 
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Model 



No Reduction 



Bias File Used 


Yes 



No 


Species 

AUCpo 

Variable 

% Cont. 

AUC P0 

Variable 

% cont. 

Paguma larvata 

0.792 

bio 22 

31 

0.794 

bio 22 

25.5 



bio_16 

20.9 


bio_16 

19.6 



bio_08 

7.2 


bio_08 

6.9 

Arctictis binturong 

0.823 

bio_02 

22.1 

0.809 

bio 22 

14 



bio 22 

17.3 


bio 02 

13.7 



bio_16 

10.6 


bio_16 

12.5 

Arctogalidia trivirgata 

0.749 

bio_02 

34.3 

0.75 

bio_02 

26.2 



bio_04 

18.1 


bio_16 

11.8 



bio_16 

13.6 


bio_04 

11.1 

Hemigalus derbyanus 

0.844 

bio_02 

24.3 

0.838 

bio 02 

19.4 



bio_04 

13.4 


bio 04 

12.9 



bio_07 

11.9 


bio_16 

12.8 

Diplogale hosei 

0.948 

bio_05 

53.1 

0.951 

bio_05 

44.3 



bio_09 

17 


bio_09 

17.8 



bio_15 

9.5 


bio_15 

12.5 

Cynogale bennettii 

0.792 

bio_04 

28.7 

0.779 

bio 04 

24 



bio_02 

21.8 


bio 02 

19.1 



bio_18 

13.3 


bio_18 

11.7 

Herpestes brachyurus 

0.844 

bio_02 

28.8 

0.836 

bio 02 

20.8 



bio 18 

17.8 


bio_18 

15.4 



bio_17 

11.3 


bio_17 

13.1 

Herpestes semitorquatus 

0.862 

bio_02 

21.1 

0.857 

bio 02 

19.1 



bio_04 

15.8 


bio_04 

14.2 



bio 22 

15.5 


bio 22 

12.8 

Neofelis diardi 

0.808 

bio_02 

30.8 

0.821 

bio_02 

22.1 



bio 12 

11.7 


bio_17 

11.5 



bio_17 

9.1 


bio 12 

10.3 

Pardofelis marmorata 

0.844 

bio_02 

28 

0.84 

bio_02 

22.9 



bio 18 

22.3 


bio_18 

20.6 



bio 22 

11.6 


bio 22 

10.5 

Catopuma badia 

0.811 

bio 22 

27.8 

0.8 

bio 22 

21.1 



bio 12 

10.7 


bio_12 

11.2 



bio_15 

8.3 


bio_15 

10.2 

Prionailurus planiceps 

0.787 

bio_04 

28.3 

0.796 

bio_04 

26.4 



bio_06 

16.6 


bio_06 

12.7 



bio_02 

8.4 


bio_03 

9.4 

Prionailurus bengalensis 

0.817 

bio_02 

35.6 

0.833 

bio 02 

23.6 



bio_04 

15.9 


bio_04 

14.2 



bio 16 

13.8 


bio 16 

14.1 
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Sensitivity analysis 2 - Spatial data filtering within 2, 5 or 10 km, respectively 


Model 



Reduction Sabah (balanced) M, 



Bias File Used 


Yes 



No 


Species 

AUCpo 

Variable 1 

% Cont. 

AUCpo 

Variable 

% cont. 

Martes flavigula 

0.735 

bio 02 

44.5 

0.73 

bio_02 

37 



bio 22 

10.8 


bio 22 

9.7 



bio 04 

8.7 


bio_04 

8.9 

Mustela nudipes 

0.806 

bio 22 

11.8 

0.806 

bio_04 

12.6 



bio 02 

11.5 


bio_19 

9.9 



bio_04 

11.3 


bio 12 

9.7 

Melogale everetti 

0.995 

bio_09 

28.7 

0.995 

bio_09 

24.8 



bio_12 

23.5 


bio 12 

24.1 



bio 17 

8.2 


bio_17 

9.2 

Mydaus javanensis 

0.912 

bio_02 

38.9 

0.908 

bio_02 

34.3 



bio_04 

32 


bio_04 

31.2 



bio_18 

9.2 


bio 18 

8.8 

Prionodon linscing 

0.781 

bio 22 

26.7 

0.782 

bio 22 

23.4 



bio_18 

22.9 


bio 18 

18.7 



bio_04 

10.8 


bio_04 

11.5 

Viverra tangalunga 

0.747 

bio 02 

25.4 

0.746 

bio_02 

22.1 



bio_04 

11.5 


bio_04 

12.7 



bio_08 

8.6 


bio_08 

7.1 

Paradoxurus hermaphroditus 

0.691 

bio_02 

30.7 

0.689 

bio_02 

27.7 



bio_04 

16 


bio_04 

16.1 



bio_18 

6.5 


bio_18 

6.7 

Paguma larvata 

0.768 

bio 22 

29.8 

0.77 

bio 22 

25.6 



bio_16 

22.8 


bio 16 

22 



bio 21 

5.9 


bio_05 

8.9 

Arctictis binturong 

0.807 

bio 22 

22.6 

0.797 

bio 22 

19.3 



bio 02 

20 


bio_02 

16.2 



bio_15 

13.6 


bio_15 

15.2 

Arctogalidia trivirgata 

0.729 

bio_02 

35.9 

0.732 

bio_02 

29.5 



bio_04 

18.2 


bio_04 

13 



bio_16 

12.4 


bio_16 

9.8 

Hemigalus derbyanus 

0.751 

bio 02 

18.7 

0.75 

bio_02 

18 



bio 04 

17.8 


bio_04 

16.1 



bio_07 

13.2 


bio_07 

10.7 

Diplogale hosei 

0.942 

bio_05 

59.6 

0.943 

bio_05 

54.2 



bio_09 

9.9 


bio_09 

9.6 



bio_15 

8 


bio_15 

9.3 

Cynogale bennettii 

0.794 

bio_04 

30.3 

0.783 

bio_04 

26 



bio 02 

20.6 


bio_02 

18.9 



bio 22 

14.1 


bio 22 

13 

Herpestes brachyurus 

0.742 

bio_02 

20.1 

0.745 

bio_02 

17.6 



bio_15 

12 


bio_15 

12.9 



bio_18 

11.2 


bio_18 

10.2 

Herpestes semitorquatus 

0.863 

bio_02 

21.4 

0.86 

bio_02 

19.7 



bio_04 

16.7 


bio_04 

15.9 



bio 22 

14.8 


bio 22 

12.3 
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Model 



Reduction Sabah (balanced) M, 



Bias File Used 


Yes 



No 


Species 

AUCpo 

Variable 1 

% Cont. 

AUCpo 

Variable 

% cont. 

Neofelis diardi 

0.772 

bio_02 

37 

0.711 

bio 02 

31.4 



bio_04 

14.3 


bio_04 

13.5 



bio 18 

7.1 


bio_18 

7.4 

Pardofelis marmorata 

0.838 

bio_02 

27.8 

0.832 

bio_02 

23.4 



bio_18 

23.3 


bio 18 

21.1 



bio 22 

11 


bio 22 

10.1 

Catopuma badia 

0.775 

bio 22 

28.9 

0.755 

bio 22 

23.8 



bio 12 

10 


bio_19 

10.9 



bio_19 

9.6 


bio 12 

10.3 

Prionailurus planiceps 

0.789 

bio_04 

25.3 

0.79 

bio_04 

23.3 



bio_06 

14.9 


bio_06 

14.1 



bio_21 

9.5 


bio_03 

10.3 

Prionailurus bengalensis 

0.701 

bio_02 

36.1 

0.703 

bio 02 

31.4 



bio_04 

21.6 


bio_04 

20.6 



bio_18 

7.5 


bio_18 

7.1 


bio_l = annual mean temperature; bio_2 = mean diurnal range (mean of monthly (max temp - min temp)); bio_3 = isothemiality (bio_2 
/ bio_7) (* 100); bio_4 = temperature seasonality (standard deviation of annual temperature *100); bio_5 = max temperature of wannest 
month; bio_6 = min temperature of coldest month; bio_7 = temperature annual range (bio_5 - bio_6); bio_8 = mean temperature of 

wettest quarter; bio_9 = mean temperature of driest quarter; bio_10 = mean temperature of warmest quarter; bio_11 = mean temperature 

of coldest quarter; bio_12 = annual precipitation; bio_13 = precipitation of wettest month; bio_14 = precipitation of driest month; bio_15 

= precipitation seasonality (coefficient of variation); bio_16 = precipitation of wettest quarter; bio_17 = precipitation of driest quarter; 

bio_18 = precipitation of warmest quarter; bio_19 = precipitation of coldest quarter. bio_21 = distance to nearest small river and water 

bodies; bio_22 = ruggedness index; bio_25 = Aspect; bio_26 = distance to medium sized rivers and water bodies; bio_27 = distance to 
main rivers and water bodies. 
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Table SM2. Spatial sensitivity analysis of species distribution models of Bornean carnivores with area calculation above fixed thresholds 
of 0.5 and 0.7 relative occurrence probability as produced by MaxEnt. Bias ‘yes’/ ‘no’ refers to inclusion of a bias file or not. ‘All’ 
variables relates to using all environmental layers as spatial input in contrast to only ‘non-correlated’ ones*. ‘No reduction’ refers to 
scenarios using all species records as input, ‘spatial filtering’ is referring to scenarios selecting only one record within a defined radius 
(M, models) and ‘balanced’ are scenarios that additionally account for the geographical bias in the input data by reducing input records 
in the state of Sabah (Mj models). 


Species 

Bias 

Variables 

Model 

%area>0.5 

%area>0.7 

Martes flavigula 

no 

all 

no reduction 

4.65 

1.54 


yes 

all 

no reduction 

13.91 

2.80 


no 

all 

Spatial filtering 

10.37 

2.60 


yes 

all 

Spatial filtering 

16.49 

3.09 


no 

all 

balanced 

18.10 

2.53 


yes 

all 

balanced 

22.82 

2.31 


no 

non-correlated 

no reduction 

6.01 

2.30 


yes 

non-correlated 

no reduction 

13.88 

3.17 


no 

non-correlated 

Spatial filtering 

12.48 

3.40 


yes 

non-correlated 

Spatial filtering 

17.10 

3.31 


no 

non-correlated 

balanced 

25.51 

1.77 


yes 

non-correlated 

balanced 

28.50 

1.57 

Mustela nudipes 

no 

all 

no reduction 

6.38 

1.44 


yes 

all 

no reduction 

14.32 

3.06 


no 

all 

Spatial filtering 

8.99 

1.98 


yes 

all 

Spatial filtering 

15.28 

2.93 


no 

all 

balanced 

8.75 

2.22 


yes 

all 

balanced 

13.95 

3.04 


no 

non-correlated 

no reduction 

8.93 

2.24 


yes 

non-correlated 

no reduction 

16.06 

3.40 


no 

non-correlated 

Spatial filtering 

10.66 

2.28 


yes 

non-correlated 

Spatial filtering 

17.15 

3.16 


no 

non-correlated 

balanced 

9.92 

2.34 


yes 

non-correlated 

balanced 

14.87 

2.95 

Melogale everetti 

no 

all 

no reduction 

0.20 

0.07 


yes 

all 

no reduction 

0.66 

0.31 


no 

all 

Spatial filtering 

0.18 

0.08 


yes 

all 

Spatial filtering 

1.15 

0.34 


no 

all 

balanced 

0.18 

0.08 


yes 

all 

balanced 

1.15 

0.34 


no 

non-correlated 

no reduction 

0.20 

0.08 


yes 

non-correlated 

no reduction 

0.75 

0.31 


no 

non-correlated 

Spatial filtering 

0.18 

0.07 


yes 

non-correlated 

Spatial filtering 

1.55 

0.45 


no 

non-correlated 

balanced 

0.18 

0.07 


yes 

non-correlated 

balanced 

1.55 

0.45 

Mydaus javanensis 

no 

all 

no reduction 

1.51 

0.48 


yes 

all 

no reduction 

4.83 

1.52 


no 

all 

Spatial filtering 

3.34 

1.36 


yes 

all 

Spatial filtering 

6.27 

2.31 


no 

all 

balanced 

6.80 

1.87 


yes 

all 

balanced 

10.04 

2.23 


no 

non-correlated 

no reduction 

2.12 

0.81 


yes 

non-correlated 

no reduction 

6.15 

2.22 


no 

non-correlated 

Spatial filtering 

4.52 

1.69 


yes 

non-correlated 

Spatial filtering 

6.94 

2.79 


no 

non-correlated 

balanced 

8.78 

1.98 


yes 

non-correlated 

balanced 

11.93 

2.22 
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Species 

Bias 

Variables 

Model 

%area>0.5 

%area>0.7 

Prionodon linsang 

no 

all 

no reduction 

9.11 

2.69 


yes 

all 

no reduction 

13.65 

3.32 


no 

all 

Spatial filtering 

9.59 

2.77 


yes 

all 

Spatial filtering 

13.80 

3.31 


no 

all 

balanced 

16.93 

2.17 


yes 

all 

balanced 

20.57 

2.04 


no 

non-correlated 

no reduction 

10.30 

3.40 


yes 

non-correlated 

no reduction 

14.07 

3.65 


no 

non-correlated 

Spatial filtering 

10.47 

3.46 


yes 

non-correlated 

Spatial filtering 

14.16 

3.68 


no 

non-correlated 

balanced 

21.62 

1.75 


yes 

non-correlated 

balanced 

24.65 

1.48 

Viverra tangalunga 

no 

all 

no reduction 

4.07 

1.32 


yes 

all 

no reduction 

14.71 

2.72 


no 

all 

Spatial filtering 

8.66 

2.46 


yes 

all 

Spatial filtering 

15.39 

2.92 


no 

all 

balanced 

20.63 

2.64 


yes 

all 

balanced 

23.55 

2.27 


no 

non-correlated 

no reduction 

5.25 

1.70 


yes 

non-correlated 

no reduction 

17.10 

2.69 


no 

non-correlated 

Spatial filtering 

10.31 

2.87 


yes 

non-correlated 

Spatial filtering 

17.02 

3.01 


no 

non-correlated 

balanced 

25.92 

1.83 


yes 

non-correlated 

balanced 

28.57 

1.43 

Paradoxurus hermciphroditus 

no 

all 

no reduction 

4.27 

1.56 


yes 

all 

no reduction 

12.29 

3.20 


no 

all 

Spatial filtering 

9.22 

2.93 


yes 

all 

Spatial filtering 

13.68 

3.32 


no 

all 

balanced 

17.41 

2.56 


yes 

all 

balanced 

22.05 

2.55 


no 

non-correlated 

no reduction 

5.36 

2.20 


yes 

non-correlated 

no reduction 

13.78 

3.31 


no 

non-correlated 

Spatial filtering 

11.60 

3.15 


yes 

non-correlated 

Spatial filtering 

15.55 

3.20 


no 

non-correlated 

balanced 

20.42 

2.58 


yes 

non-correlated 

balanced 

24.36 

2.02 

Paguma larvata 

no 

all 

no reduction 

9.69 

2.03 


yes 

all 

no reduction 

14.92 

3.34 


no 

all 

Spatial filtering 

10.42 

2.22 


yes 

all 

Spatial filtering 

14.96 

3.14 


no 

all 

balanced 

11.84 

3.08 


yes 

all 

balanced 

15.88 

3.17 


no 

non-correlated 

no reduction 

13.20 

2.71 


yes 

non-correlated 

no reduction 

18.53 

3.62 


no 

non-correlated 

Spatial filtering 

14.08 

3.21 


yes 

non-correlated 

Spatial filtering 

18.65 

3.58 


no 

non-correlated 

balanced 

13.10 

3.79 


yes 

non-correlated 

balanced 

18.18 

3.73 
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Species 

Bias 

Variables 

Model 

%area>0.5 

%area>0.7 

Arctictis binturong 

no 

all 

no reduction 

4.58 

1.12 


yes 

all 

no reduction 

13.16 

2.71 


no 

all 

Spatial filtering 

7.76 

1.98 


yes 

all 

Spatial filtering 

14.60 

2.80 


no 

all 

balanced 

12.19 

2.99 


yes 

all 

balanced 

16.01 

3.36 


no 

non-correlated 

no reduction 

6.32 

1.58 


yes 

non-correlated 

no reduction 

15.37 

2.72 


no 

non-correlated 

Spatial filtering 

10.33 

2.59 


yes 

non-correlated 

Spatial filtering 

15.10 

2.90 


no 

non-correlated 

balanced 

13.30 

3.68 


yes 

non-correlated 

balanced 

16.04 

3.75 

Arctogalidia trivirgata 

no 

all 

no reduction 

10.11 

2.38 


yes 

all 

no reduction 

17.24 

3.36 


no 

all 

Spatial filtering 

11.02 

2.25 


yes 

all 

Spatial filtering 

17.44 

3.24 


no 

all 

balanced 

13.38 

1.80 


yes 

all 

balanced 

20.30 

2.27 


no 

non-correlated 

no reduction 

11.97 

2.69 


yes 

non-correlated 

no reduction 

18.46 

3.29 


no 

non-correlated 

Spatial filtering 

12.86 

2.68 


yes 

non-correlated 

Spatial filtering 

18.91 

3.22 


no 

non-correlated 

balanced 

15.53 

2.07 


yes 

non-correlated 

balanced 

21.80 

2.29 

Hemigalus derbyanus 

no 

all 

no reduction 

5.01 

1.55 


yes 

all 

no reduction 

14.37 

2.77 


no 

all 

Spatial filtering 

8.96 

2.66 


yes 

all 

Spatial filtering 

14.80 

3.02 


no 

all 

balanced 

17.77 

2.37 


yes 

all 

balanced 

23.60 

2.20 


no 

non-correlated 

no reduction 

6.53 

2.02 


yes 

non-correlated 

no reduction 

14.50 

2.84 


no 

non-correlated 

Spatial filtering 

11.60 

2.98 


yes 

non-correlated 

Spatial filtering 

16.24 

3.02 


no 

non-correlated 

balanced 

20.01 

2.88 


yes 

non-correlated 

balanced 

25.18 

2.44 

Diplogale hosei 

no 

all 

no reduction 

2.50 

0.51 


yes 

all 

no reduction 

4.38 

0.94 


no 

all 

Spatial filtering 

2.81 

0.58 


yes 

all 

Spatial filtering 

4.32 

0.91 


no 

all 

balanced 

1.87 

0.46 


yes 

all 

balanced 

3.48 

0.84 


no 

non-correlated 

no reduction 

3.08 

0.61 


yes 

non-correlated 

no reduction 

4.84 

0.92 


no 

non-correlated 

Spatial filtering 

3.42 

0.71 


yes 

non-correlated 

Spatial filtering 

4.80 

0.97 


no 

non-correlated 

balanced 

2.36 

0.51 


yes 

non-correlated 

balanced 

3.92 

0.81 
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Species 

Bias 

Variables 

Model 

%area>0.5 

%area>0.7 

Cynogale bennettii 

no 

all 

no reduction 

2.49 

1.21 


yes 

all 

no reduction 

6.34 

2.34 


no 

all 

Spatial filtering 

2.58 

1.24 


yes 

all 

Spatial filtering 

5.90 

2.27 


no 

all 

balanced 

15.70 

2.25 


yes 

all 

balanced 

18.21 

2.27 


no 

non-correlated 

no reduction 

3.33 

1.50 


yes 

non-correlated 

no reduction 

6.39 

2.65 


no 

non-correlated 

Spatial filtering 

3.21 

1.50 


yes 

non-correlated 

Spatial filtering 

6.05 

2.47 


no 

non-correlated 

balanced 

22.77 

1.67 


yes 

non-correlated 

balanced 

24.89 

1.61 

Herpestes brachyurus 

no 

all 

no reduction 

3.50 

1.19 


yes 

all 

no reduction 

11.17 

2.56 


no 

all 

Spatial filtering 

11.91 

2.69 


yes 

all 

Spatial filtering 

16.60 

2.82 


no 

all 

balanced 

16.49 

2.13 


yes 

all 

balanced 

20.13 

2.18 


no 

non-correlated 

no reduction 

4.09 

1.58 


yes 

non-correlated 

no reduction 

13.03 

2.72 


no 

non-correlated 

Spatial filtering 

14.57 

3.40 


yes 

non-correlated 

Spatial filtering 

18.90 

3.55 


no 

non-correlated 

balanced 

18.97 

2.12 


yes 

non-correlated 

balanced 

21.92 

2.26 

Herpestes semitorquatus 

no 

all 

no reduction 

2.24 

1.08 


yes 

all 

no reduction 

7.67 

2.10 


no 

all 

Spatial filtering 

2.51 

1.21 


yes 

all 

Spatial filtering 

7.33 

2.18 


no 

all 

balanced 

11.54 

1.85 


yes 

all 

balanced 

14.05 

2.00 


no 

non-correlated 

no reduction 

2.97 

1.24 


yes 

non-correlated 

no reduction 

8.32 

2.41 


no 

non-correlated 

Spatial filtering 

3.44 

1.45 


yes 

non-correlated 

Spatial filtering 

7.73 

2.45 


no 

non-correlated 

balanced 

12.84 

1.54 


yes 

non-correlated 

balanced 

15.94 

1.64 

Neofelis diardi 

no 

all 

no reduction 

3.16 

1.20 


yes 

all 

no reduction 

12.28 

2.59 


no 

all 

Spatial filtering 

6.85 

2.02 


yes 

all 

Spatial filtering 

13.30 

2.95 


no 

all 

balanced 

19.92 

2.39 


yes 

all 

balanced 

24.41 

2.09 


no 

non-correlated 

no reduction 

4.80 

1.67 


yes 

non-correlated 

no reduction 

14.37 

2.85 


no 

non-correlated 

Spatial filtering 

8.75 

2.43 


yes 

non-correlated 

Spatial filtering 

15.35 

2.97 


no 

non-correlated 

balanced 

25.32 

2.01 


yes 

non-correlated 

balanced 

28.44 

1.80 
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Species 

Bias 

Variables 

Model 

%area>0.5 

%area>0.7 

Pardofelis marmorata 

no 

all 

no reduction 

3.92 

1.52 


yes 

all 

no reduction 

9.36 

2.98 


no 

all 

Spatial filtering 

4.15 

1.64 


yes 

all 

Spatial filtering 

8.87 

2.92 


no 

all 

balanced 

13.14 

2.33 


yes 

all 

balanced 

17.71 

2.47 


no 

non-correlated 

no reduction 

5.11 

1.98 


yes 

non-correlated 

no reduction 

10.09 

3.03 


no 

non-correlated 

Spatial filtering 

5.44 

2.11 


yes 

non-correlated 

Spatial filtering 

9.70 

3.00 


no 

non-correlated 

balanced 

14.24 

2.46 


yes 

non-correlated 

balanced 

19.01 

2.48 

Catopumci badia 

no 

all 

no reduction 

3.53 

1.07 


yes 

all 

no reduction 

9.87 

2.14 


no 

all 

Spatial filtering 

4.17 

1.28 


yes 

all 

Spatial filtering 

9.47 

2.15 


no 

all 

balanced 

15.74 

1.81 


yes 

all 

balanced 

17.93 

1.84 


no 

non-correlated 

no reduction 

5.71 

1.73 


yes 

non-correlated 

no reduction 

12.66 

2.24 


no 

non-correlated 

Spatial filtering 

7.49 

1.94 


yes 

non-correlated 

Spatial filtering 

12.63 

2.34 


no 

non-correlated 

balanced 

21.85 

1.81 


yes 

non-correlated 

balanced 

24.66 

1.78 

Prionailurus planiceps 

no 

all 

no reduction 

6.67 

2.03 


yes 

all 

no reduction 

13.24 

3.25 


no 

all 

Spatial filtering 

7.51 

2.22 


yes 

all 

Spatial filtering 

13.03 

3.10 


no 

all 

balanced 

10.70 

2.26 


yes 

all 

balanced 

15.44 

2.63 


no 

non-correlated 

no reduction 

10.37 

2.94 


yes 

non-correlated 

no reduction 

16.05 

3.70 


no 

non-correlated 

Spatial filtering 

11.61 

3.15 


yes 

non-correlated 

Spatial filtering 

15.95 

3.59 


no 

non-correlated 

balanced 

13.92 

2.18 


yes 

non-correlated 

balanced 

17.38 

2.40 

Prionailurus bengalensis 

no 

all 

no reduction 

5.09 

1.63 


yes 

all 

no reduction 

15.74 

2.81 


no 

all 

Spatial filtering 

9.65 

2.70 


yes 

all 

Spatial filtering 

15.62 

3.01 


no 

all 

balanced 

16.51 

2.13 


yes 

all 

balanced 

21.89 

2.07 


no 

non-correlated 

no reduction 

7.31 

2.35 


yes 

non-correlated 

no reduction 

16.14 

3.51 


no 

non-correlated 

Spatial filtering 

10.70 

3.09 


yes 

non-correlated 

Spatial filtering 

16.54 

3.16 


no 

non-correlated 

balanced 

20.82 

1.77 


yes 

non-correlated 

balanced 

25.86 

1.71 


*Non-correlated variables comprise: bio_l, 4, 7, 13, 14, 15, 22, 26 and 27 (for description see table subtext SMI; for correlation matrix: 
see Supplementary Material in Kramer-Schadt et al., 2013). 
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